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Abstract. We prove that the distribution functions of magnetization and spin glass 
order parameter coincide on the Nishimori hnc in the phase diagram of the ± J Ising 
model in any dimension. This implies absence of replica symmetry breaking because the 
distribution function of magnetization consists only of two delta functions, suggesting 
the same simple structure for the distribution of spin glass order parameter. It then 
follows that the mixed (glassy) phase, where the ferromagnetic order coexists with 
complex phase space, should lie, if any, below the Nishimori line. We also argue that 
the AT line to mark the onset of RSB with a continous distribution of the spin glass 
order parameter, if any again, would start with an infinite slope from the multicritical 
point where paramagnetic, ferromagnetic and spin glass phases merge. 



INTRODUCTION 

Existence and characteristics of the spin glass phase are actively investigated for 
finite-dimensional random spin systems. Closely related is the problem of the mixed 
(glassy) ferromagnetic phase. If the mean-field picture applies to finite-dimensional 
systems, the ferromagnetic phase would split up into two regions, one with a simple 
structure (the replica-symmetric (RS) phase in the mean-field framework) and the 
other with a complex phase space (replica-symmetry broken (RSB) state). The 
boundary between these two phases would be the finite-dimensional counterpart of 
the AT (de Almeida-Thouless) line found for the infinite-range SK (Sherrington- 
Kirkpatrick) model. 

Investigations of these and related problems are almost exclusively carried out 
currently by numerical methods because of lack of reliable analytical methods. We 
show in the present contribution that simple symmetry arguments lead to a strong 
constraint on the possible location and shape of the AT line in the phase diagram 
of the Ising spin glass on an arbitrary lattice with arbitrary range of interactions. 



More precisely, it is possible to prove that there is nothing hke RSB on the Nishimori 
hne in the phase diagram and that the AT line, if it ever exists, should start as 
a vertical line from the multicritical point where paramagnetic, ferromagnetic and 
spin glass phases merge. Our method is an application of the gauge theory which 
has been used to derive exact energy and many other exact/rigorous results. 

GAUGE THEORY OF THE ISING SPIN GLASS 

Let us consider the ±J Ising model {Si — ±1) with the Hamiltonian 

H ^ -J^JijSiSj, (1) 

(ij) 

where the sum is over pairs of sites on an arbitrary lattice with arbitrary range of 
interactions. The exchange interactions are quenched random variables with the 
distribution function 

P( J,,) = p6{J^j - J) + (1 - p)6{J,, + J). (2) 

The Hamiltonian (1) is invariant under gauge transformation 

Jij ^ Jijf^i^^ji Si > Sidi, (3) 

where cij is an arbitrarily fixed Ising spin. 

The gauge invariance of the Hamiltonian leads to a number of exact and/or 
rigorous results [1,2]. For example, consider the internal energy 

E{K,K,) = [{H)k]k^, (4) 

where the inner triangular brackets denote a thermodynamic average and the outer 
square brackets represent the configurational average at a fixed value of p. The 
suffixes denote the temperature and p through K = (3 J = J/ksT and Kp = 
I log]9/(l — p). The exact value of the internal energy can be calculated explicitly 
under the constraint Kp — K to give 

E{K,K) = -NeJianhK, (5) 

where Nb is the number of bonds on the lattice under consideration. The constraint 
K = Kp defines a curve in the phase diagram shown dashed in Figure 1, which is 
often called the Nishimori line. 

There are many other results including an upper bound on the specific heat 

kBT'CiK,K)<—^, (6) 
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FIGURE 1. Phase diagram of the ±J model and the Nishimori Une (dashed). 



[{SiSj)K]K 

and a correlation inequality 



{SiSj)\ 



{SiSj 



(7) 



(8) 



The correlation identity (7) means that the Nishimori line {Kp = K) avoids the spin 
glass phase because the left-hand side reduces to the square of the ferromagnetic 

long range order as the distance between i and j tends to infinity whereas 
the right-hand side approaches the square of the spin glass order parameter in 
the same limit. The spin glass phase is characterized by g > 0,m = which is 
inconsistent with the relation m — ±g. The correlation inequality (8) places a 
constraint on the possible shape of the boundary between the ferromagnetic and 
spin glass phases: This boundary can either be vertical or re-entrant in the phase 
diagram. The ferromagnetic phase is not allowed to lie below the spin glass phase 
[1,2]. 

DISTRIBUTION OF ORDER PARAMETERS 

The gauge theory summarized above applies to derive an identity between the 
distribution functions of the ferromagnetic and spin glass order parameters on the 
Nishimori line [2,3]. 



Definitions and the main result 

The distribution function of magnetization is defined as 



Pm{x; K) 



(9) 



To define the distribution function of the spin glass order parameter, it is convenient 
to introduce two replicas of the same system with spins {Si\ and {(7^} and the 
couplings Ki and K21 respectively: 




FIGURE 2. Distribution of magnetization (left) and the SG order parameter (right). 



-PiH{S)p~l32H{a) 



(10) 



where Ki = f3iJ and K2 = (32 J ■ The distribution of the spin glass order param- 
eter is defined in terms of two rephcas with the same temperature, Pq{x;K) = 
Pq{x; K, K). Note that the distributions Pm{x; K), Pq{x; K, Kp) and Pg{x;K) are 
all functions of p as well. 

It is well established that the distribution of magnetization has a simple structure 
with two delta functions as depicted in the left part of Figure 2. The spin glass order 
parameter, on the other hand, has a more complex structure. For example, if there 
exists a full RSB, the distribution has a continuous part as depicted schematically 
in the right part of Figure 2. Or, it may have delta peaks at a few additional 
locations in the case of finite-step RSB. In any case, the complexity of the system 
is closely related with the structure of the distribution function Pq{x) different from 
that of Pm{x). 

Our main result is the following relation: 

Pm{x;K) = Pq{x;K,Kp) (11) 

which holds for any combination of x, K and Kp. In particular, if wc set Kp = K. 
Equation (11) reduces to the relation between the usual distribution functions of 
magnetization and spin glass order parameter: 

Pr,,{x;Kp)^Pq{x;Kp). (12) 

This immediately proves that the distribution function of the spin glass order pa- 
rameter on the right-hand side has a simple structure represented by the magneti- 
zation distribution on the left-hand side ii K = Kp. 

The simplicity of the exact energy (5) and a few other results [2,1] suggest that 
the phase space of the system would not be complicated. The above result is 
regarded as definite evidence to support this anticipation. 

It also follows from differentiation of Equation (11) that the derivatives of the 
distribution functions Pm{x'i K) and Pq{x] K) with respect to K and p vanish when 
K^Kp. 



FIGURE 3. Phase diagrams compatible (left) and incompatible (right) with the relations proved 
in the text when the mixed phase is characterized by a continuous part in Pq{x). 



— P,(x;X) = 0, -P,{x;K)^0 (13) 

for almost all x when K = Kp. These relations have been derived from the fact 
that the distribution of magnetization Pm{x] K) consists just of two delta functions 
and therefore its derivatives with respect to the parameters K and p vanish for 
almost all x. Equation (13) then follows from (11). 

The relation (13) shows that the distribution function of the spin glass order 
parameter remains vanishing at almost all x if the parameters K and p deviate 
infinitesimally from their values satisfying K = Kp. Thus the structure of Pq{x] K) 
does not develop a continous part under an infinitesimal deviation from the Nishi- 
mori line. An interesting consequence is that the mixed phase with a continuous 
part in Pq{x)^ if it starts to exist at the multicritical point, does not have a finite 
extension under a very small deviation from the multicritical point. Stated other- 
wise, the phase boundaries (the AT line and the boundary between the mixed and 
spin glass phases) should merge smoothly as they approach the multicritical point 
as depicted on the left part of Fig. 3. A natural generic situation would be that 
both of these boundaries are vertical around the multicritical point as in the case 
of the SK model. It should be noted that this argument does not apply when the 
RSB in the mixed phase is characterized by several delta functions, in which case 
the boundary with the structure like the one on the right part of Fig. 3 does not 
violate (13). 



Outline of the proof 

The proof of the relation (11) is not very difficult [3]. We start from writing out 
the configurational average of the distribution function of magnetization 

^"^^"^'""^ (2coshi^,)^B^^^;e;i/ ^^e^E-.^«^. ■ ^'^> 



The weight of configurational average (p or 1 — p for each bond) is taken care of by 
the factor e^^S"^*^ [1,2]. The gauge transformation (3) changes the above equation 
into the form 

^m{X,l^) (2coshXp)^s ^^^^tii} ^^e^E-*.-5i5, 

(15) 

It is relatively straightforward to see that the right-hand side is equal to 
Pg{x;Kp,K) by summing it up over all possible {o-j}, dividing the result by 2^, 
and finally inserting an identity 1 = J2a e^p^'^^^'^''^^ / J2a e^^^'^*^'^*'^^ after the sum- 
mation symbol over {Tij}. 



CONCLUSION 



We have proved that the AT line, if any, should lie below the Nishimori line in 
the phase diagram of the ±J Ising model on an arbitrary lattice. It has also been 
argued that the mixed phase with a continous part in the distribution function of 
the spin glass order parameter, if any again, does not have a finite extension around 
the multicritical point. These results have been derived by a simple application of 
gauge transformation. 
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